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Recent lattice studies exhibit infrared finite effective QCD charges. Corresponding gluon prop- 
agator in Landau gauge is finite and nonzero, suggesting a mechanism of dynamical gluon mass 
generation is in the operation. In this paper, the analytical continuation of the Euclidean (space- 
like) Pinch Technique-Background Field Method (PT-BFM) solution of Schwinger-Dyson equation 
for gluon propagator to the timelike region of q 2 is found. We found the continuation numerically 
showing good agreement with a generalized Lehman representation for small Schwinger coupling. 
The associate non-positive spectral function has an unexpected behavior. Albeit infrared Euclidean 
space solution naively suggests like single scale "massive" propagator, the obtained spectrum of 
gluon propagator does not correspond to the delta function at single scale q — m , instead more 
f\| ' possible singularities are generated. The pattern depends on the details of assumed Schwinger 

mechanism: for stronger coupling there are few maxima and minima which appear at the scale A, 
while for perturbatively small Schwinger coupling the spectral function shows up two narrow peaks: 
04 ' particle and ghost excitation, which have mutually opposite signs. 
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I. INTRODUCTION 

Infrared behavior of Green's functions (GFs) of pure Yang-Mills theories has been intensively studied in the last 
decade. Even though GFs are gauge fixing dependent objects and thus they do not represent physical observables, it 
is believed they include reliable nonperturbative information about confinement (and chiral symmetry breaking when 
quarks are included). Various solutions for GFs have been obtained by solving Schwinger-Dyson equations (SDEs) and 
by lattice calculations. Recently, study of SDEs offer two scenarios which are distinguished by the infrared behavior 
i— i' of GFs. The so called scaling solution exhibits power law momentum behavior with an infrared exponents [H-Q which 
lead to infrared vanishing Landau gauge gluon propagator A(0) = and correspondingly infrared enhanced ghost 
propagator. For a SDEs treatment in QCD let us refer (HQ. 
ly-j , On the other side, recent lattice calculations [fj|-[l8j in conventional Landau gauge support the so called decoupling 
• solution. In this scenario the gluon propagator A(0) ^ 0, but stays finite. Also, it is notable, the ghost propagator 
is not enhanced in the infrared and thus remains semi-perturbative. Such decoupling solutions, called sometimes 

■ massive solutions were already proposed in [19|-|21| and being based on attractive theoretical frameworks of Pinch 
; Technique and Background Field Method (PT-BFM) they attract many new attentions these days [H-IH, [UrS 51 ■ 

Comparing to standard gauge fixed set of GFs it has been proved that using Pinch Technique- one can rearrange 
t-H ' usual gauge fixed GFs (in any gauge) in a way that they do not depend on the gauge fixing parameters and PT GFs 
satisfy original WTI an not more complicates Slavnov- Taylor identities, for the topical review see [27|. While scaling 
solution would be compatible with Gribov-Zwanziger confinement scenario, it is believed that the decoupling solutions 
is related to a dynamical generated gluon mass m? ~ A _1 (0), however confinement mechanism of associated massive 

■ gluons remains to be determined. 

The truncation of such reorganized SDEs with WTI satisfying internal vertices is called "PT-BFM inspired" or 
simply PT-BFM gluon SDE. As a starting point we use the model originally solved in Euclidean space and for certain 
set of free parameters presented in the paper [22j • The dynamical mass generation is a BFM generalization of the well 
known Schwinger mechanism, hereby singular vertex appears in quantum loops only and not due to the Goldstone 
modes, which do no take place here as the gauge symmetry is exactly preserved. However as we do not solve more 
complicated SDE for the gluon vertices, we freely accommodate modeled transversal part in its original form [22j . 
Having solved the problem in the Euclidean space we formally perform analytical continuation of the equation to the 
timelike region, i.e. to the Minkowski space, where it should provide analytical continuation as a solution. 

In BFM scheme, the expected Schwinger mechanism gives rise to the infrared gluon mass which naively could not 
be so far from the value of A _1 / 2 (0), while the true physical particle mass is determined as the pole mass of the full 
propagator. For a stable unconfined particle it can be extracted as a real solution of equation: A _1 (gr 2 = m p (q 2 )) = 0. 
Gluon is experimentally unobservable as a particle and thus very likely free plane wave solution does not exists. If 
there is no real pole in the gluon propagator then a propagation of on-shell gluon is indeed impossible. What is 
the singularity structure of the propagator of expectingly confined excitations has not been known. Instead of this, 
some attempts to determine (not completely understood) gluon mass scale were performed on various theoretical and 
phenomenological basis |29l433 | giving us rough and simplified estimate m g = Aqcd- 
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In this paper we study Minkowski space continuation of already considered PT-BFM solution 22j ,which although 
well defined in the entire Minkowski space, was numerically performed in Euclidean space only. In reality, we did 
not get just one solution of continued PT-BFM gluon SDE but many, among them we have to choose the correct 
solution. For this purpose it is fully legitimate to check the analytical properties which are equivalent to the analyticity 
Stieltjes transformable function (generalized Lehman representation). This is fully consistent treatment as the Lehman 
representation was explicitly used during the derivation of the original Euclidean space SDE. Due to this reason we will 
abandon all the solutions which largely do not fit dispersion relation dictated by Lehman representation. Optimizing 
the analyticity as far possible we got a few scenarios, albeit they are quite sensitive to the details of the model. 
Especially, and quite interestingly it allows to answer what could be a physical mass of the gluon in the pure Yang- 
Mills theory Neither we found single solution A~ 1 (q 2 = m 2 ) = for some pole mass m p . Actually, solutions we have 
found lie between two extreme cases: first is the solution with two particle like singularities with opposite signs and 
the second case, where there is no solution for the physical mass shell. In later case the associated spectral function 
is an oscillatory one with a few number of relatively smooth minima and maxima. The later case we can interpret 
as a spectrum of confined gluon, while the first kind of solutions is very likely an artifact of approximation made. In 
both cases the solution suggests the first branch point is located in the origin of complex q 2 plane. 

Our convention for Minkowski metric tensor reads: — diag(+l,—l,—l,—l). Minkowski momenta are not 
labeled in our notation, while we always write E when we specify the Euclidean momenta, i.e. for instance q 2 E = —q 2 
for some spacelike momentum q. We use simple letter d for the PT-BFM gluon propagator in Minkowski and c?£ in 
Euclidean space, while symbol A is used for a conventional propagator defined in given gauge. In the next Sections II 
and III we describe necessary ingredients of the original PT-BFM SDE and perform formal " analytical" continuation. 
In Section IV we continue with RG improved SDEs for which we present numerical solutions. We discuss what 
happens to the solutions in various cases, we also discuss necessary amount of numerics, which is the key tool to get 
a correct and stable results in Minkowski space. In the Section V we mention some outcomes for continuation of the 
lattice gluon propagator and we conclude and discuss some open questions in Section VI. 



II. ANALYTICAL CONTINUATION I, PT-BFM GLUON SDE 



In quantized gauge theory one needs to fix gauge in order to be able to calculate GFs. Usual approach is Fadeev- 
Popov [33| gauge fixing procedure which makes perturbation calculation particularly feasible, however GFs become 
gauge fixing dependent. Pinch Technique is the S-matrix based construction of gauge fixing independent GFs. After 
reorganization of GFs (calculated at any gauge) one get new GFs which satisfy Ward Identities instead of usual more 
complicated Slavnov- Taylor identities. It has been proved to all order of the perturbation that Pinch Technique is 
equivalent to the theory quantized by Background Field Method at Feynman gauge. Guiding by the principles of 
Pinch Technique, the PT-BFM gluon SDE has been recently studied in [22l. I29I] . 

To get infrared finite gluon propagator (and running coupling), the gluon propagator must loose its perturbative 
1/q 2 pole through the Schwinger mechanism in Yang-Mills theory [43]. It underlies on the assumption of infrared 
singular gluon vertex in a way it leaves polarization tensor transverse (gauge invariant). In this paper we simply use 
the derived PT-BFM equation in [22j], where Schwinger mechanism is employed through the simple Ansatz for the 
improper (two leg dressed) three gluon PT-BFM dressed vertex 
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where T^ a/3 satisfies tree level WTI and d is scalar function related to the all order PT-BFM gluon propagator which 
in Landau gauge reads 
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and satisfies generalized Lehman representation 
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and is the rest of the three gluon improper vertex which is not specified by gauge invariance. The essential feature 
of the vertex T ua p is that apart the structure dictated by WTI it also includes 1/q 2 pole term which gives rise to 
infrared finite solution. For this purpose the following form 
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has been proposed in [22j ]. This vertex is transverse in respect to q (q.T = 0) and it respects Bose symmetry to two 
quantum legs interchange as its origin is due to the quantum loops. In this respect there are no singularities associated 
with external gluons, avoiding thus singular footprint in (however pcrturbative) any S-matrix element. 
Finally, after the renormalization, it leads to the following form of linearized SDE in Euclidean space: 
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where the second line arises due to the Ansatz for the gluon vertex (j2.4[) and K is the renormalization constant. Thus 
the strength of the dynamical mass generation is triggered through the adopted coupling constants c\ , C2 , C3 which is 
fully equivalent to the introduction of (in principle arbitrary) constant 7 and infrared value eT^^O) (for completeness 
recall that in the paper [22| d(0) has be calculated since one of a was fixed by hand). 

The analytical assumption (|2.3[) is the right key for finding a correct continuation in Minkowski space. All the 
singularity structure is encoded in the distribution p. Thus if one knows somehow the essential type of singularity 
(the example is simple, i.e. p — 5{x — a) in this case) then one can estimate the positions of induced singularity 
structure, e.g. the position of the first branch points (thresholds are examples generated by simple poles of internal 
propagators). Consequently , after making formal analytical continuation of SDE one can evaluate the imaginary and 
the real part of propagator function at domain where the function d is continuous. This procedure actually works 
numerically for the models, where the continuous function part of p is smooth enough and associated principal value 
integrals can be performed with high numerical accuracy, see [3~i| for a review. However when the coupling constant 
is large enough, one is faced to an oscillatory behavior of the spectral function p [35j-|37[ and the SDE as an integral 
equation for p becomes unstable. We found the system of such equations extremely unstable in the case of PT-BFM 
gluon SDE as well and due to this reason we develop different strategy for this purpose. 

Due to a simple linear structure of the SDE (12.51) the analytical continuation is very straightforward. Recall the 
definition here 

d E {q 2 E ) = -d(-q 2 ) for q 2 < , (2.6) 

where now d is Minkowski space PT-BFM gluon propagator. Putting q 2 E — > —q 2 and considering positive (timelike) 
q 2 continuation we get the SDE for continued gluon propagator: 
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where the substitution z — > — z was made and also the both sides were multiplied by —1. 

While the Eq. (12.71) is formally identical to its spacelike counterpartner , the main difference should be stressed: 
The function d on the both sides is complex function, the one on the rhs. of SDE is integrated over the region where 
all singularities are located. Meaning of Feynman it is as usual, it is responsible for the generation of absorptive part 
of d at the timelike region, however here, we keep it small and finite in our numerical treatment. Having d complex, 
the equation (|2.7p represents two coupled real equations for an imaginary and the real parts of d. Thus the presence 
of small imaginary is crucial for obtaining a correct Minkowski solution. Clearly switching off imaginary factor we 
would repeat Euclidean space solution again, albeit now for timelike q 2 . Clearly such solution do not complete the 
full solution in Minkowski space since the timelike solution of Eq. (|2.7p is purely real and contradicts the assumption 
made. Also note, such solution has an infrared discontinuity d~ 1 (0 + ) — d _1 (0~) = 2A 2 (left limit is the Euclidean) 
, while for a continuous solution we need to adjust d (0+) = —A 2 . At this place, we can mention we have found a 
real solutions of Eq. (12. 7p for d~ 1 (0 + ) < , which includes the case d _1 (0 + ) = —A 2 as well. Up to the later case all 
of them posses second order discontinuity, while for the choice d~ 1 (0 + ) = —A 2 we also got the real solution but with 
the first order singularity only. Of course, these real solutions cannot be interpreted as the analytical continuation 
which we are looking for since the reality of d completely contradicts our assumptions. Much interesting for us is the 
following observation: for mentioned solutions we get the second and the first order singularity respectively, which is 
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located at q 2 = 0. It is the motivation to look, if not yet an evidence, for the complex solution with the branch point 
located at the origin of q 2 complex plane. 

We argue here that the solution of Eq. (|2.7[) is not unique, there must be more solutions in addition, at least one 
is complex and respects Lehman representation and corresponds thus to the continuation to the Minkowski timelike 
domain of q 2 , this later has the branch point located at the origin of Minkowski q 2 = 0. There exist probably 
more (infinite number is not excluded) solutions in other Riemann sheets if d is a multivaluable analytical functions. 
The required solution is represented by the function which starts to be complex from the origin in accordance with 
the introduction of trivial boundary in the integral expression (|2.3[) . Actually we got numerical evidence for such 
statement by our numerical search. 

Before discussing more details, we improve the linearized SDE in a way it gives correct UV asymptotic behavior, 
for which case we also present the numerical solution. 

III. RENORMALIZATION GROUP IMPROVED GLUON SDE 

Solving the SDE perturbatively by taking d(z) —tl/z one get the correct 1 loop perturbative limit 

d{x) ^-* x{l + bg*ln{ X /f) - (3 - 1} 

However, when solving SDE selfconsistently then the solution of (|2.7p leads to the slower decrement with growing 
momenta q 2 , instead of log suppression we would get (1 + bg 2 ln(x/ n 2 )) 1 / 2 (see [22[ for the details) behavior of the 
gluonic form factor, which is believed is the prize of the gauge technique and associated renormalization simplifications. 
In order to restore the correct RG behavior of the SDE solution one needs to replace 

g 2 d{z) -^—Mz) (3.2) 

9uv( z ) 

in the integrands of SDE, where g 2 rv (z) is some function with the perturbation theory running coupling ultraviolet 
asymptotic. 

To accomplish this we take the function gjjv in the form 

g ^ (z) = Mn(e+£) - (3 ' 3) 

Up to an required asymptotic it is smooth function enough, which ensures unwanted infrared modification of the SDE 
kernel. 

As a first necessary step we obtain the solution of PT-BFM SDE in the spacelike region and we solve the SDE 
in the Euclidean space. Incorporating substitution (|3.2I) . fitting the renormalization constant K = 1 we solved the 
following equation 
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As in the previous case of Eq. (12.71) , two parameters stay completely free in given truncation. For simplicity and 
quite naturally we take the infrared gluon mass identical with AqCD. 

d- 1 (0) = A 2 . (3.5) 

and then looking what happens when 7 is varied. We have found there is little dependence of the solution even when 
we change 7 about order of magnitude, for the solutions see Fig. [TJ 

After a slight rearrangement, the RG improved PT-BFM SDE (|3.4[) has been solved by the iteration. The iteration 
procedure shows up very fast convergence, providing 20-30 iteration steps are enough to achieve vanishing difference 
between the two followed iterations. Since the Euclidean correlator is quite smooth function thus the resulting 
solutions are numerically stable. Thanks to the I-dimensionality of the problem we can use a grid with large number 
of mesh points giving easily 3 to 6 digit accuracy without a special treatment of the upper boundaries for which we 
use simple step function. In our treatment we do not define a new gluon form factor neither gluon dynamical mass as 
this cannot be done without ambiguity. In figure [T] we plot the resulting solution for the dressing function F = q 2 d. 
The comparison with the I-loop perturbative running coupling asymptotic is shown in figure [3J 
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IV. ANALYTICAL CONTINUATION II, RESULTS FOR RENORMGROUP IMPROVED PT-BFM SDE 

If gluon were a stable massive particle we could get delta distribution S(lu — p 2 ) as a first substantial singularity in 
gluon spectral function. Further there could be a second singularity at 4p 2 where the continuous part of spectra usually 
starts. As the massless ghosts were omitted we could get further singularities at 9/z 2 , 16/i 2 , ... all these singularities 
exhibit ourself as a finite cusps in the real part of d, since p has no derivatives at these points and all of them could 
be located on the cut in complex plane which is identical with the real positive semiaxis of Minkowski q 2 . 

However, instead of dealing with stable vector particle, we are dealing with confined gluon field. Let us imagine 
a process when "high energy" gluon is emitted with some virtuality q 2 . It starts to loose its energy through the 
pair creation and bremsstrahlung of soft gluons, all these colored intermediate states are finally (and simultaneously) 
neutralized in a color singlet glueballs (hadrons, when real QCD with quarks is considered). Such color neutralization 
can happen only through the interaction with other colored partons- gluons and quarks. Gluon never escapes as a 
free on shell particle, no matter what the typical gluon mass scale would be. From this picture it is apparent: the 
timelike structure can be quite complicated as many "intermediate states" contributes significantly because of strong 
coupling. The emission and absorption of quanta of the fields are related with threshold singularities when unconfincd 
particles are considered. Here, for the case of gluons, we expect something similar, but perhaps qualitatively different: 
the absorption and emission of confined object could reflect their short life t ~ 1 /Aqcd , we expect this reflection is 
encoded in the singularity structure typical for confinement. Thus considering spectral function one can really expect 
kind of cusps instead of delta function and some more smooth non-monotonic behavior instead of threshold cusps. 
Constrained by construction here, we do not expect their positions in a complex plane away from a positive real axis. 

In any case, assuming spectral function has a continuous part we must get standard dispersion relation for 



where P. stand for principal value integration. 

Without knowledge of branch points and singularities in p one is faced with problem of performing limiting ana- 
lytical continuation to the borderline of analyticity- to the cut, where at (hopefully) isolated points we assume d has 
aforementioned singularities of not yet specified type. Not necessarily but likely, they can be related with a branch 
points of multivaluable analytical function and very likely the number of them increases with growing q 2 . It is our 
assumption that all the information about singularities is correctly captured in Minkowski space continued SDE for 
gluon. Of course , one can expect that various approximation, which are necessary complements of given truncation 
of SDEs system has an impact on the solution, e.g. on the position, number of branch points and even the analyticity 
can be destroyed. Finally remind the well known examples: this is the Dirac delta function in spectrum which im- 
plies pole in the propagator, similarly Heaviside step function 6(x — a) gives log(q 2 — a) divergence and for instance 
9(1 — T /x)yf(\ — T/x) gives usual particle production threshold at T. That it is easy to write down a plethora of 
arbitrarily moderated singularities, which are not familiarly known from particle physics is quite obvious. 

In previous discussion we offer some arguments which are crucial for the method of the solution. Let us summarize 
and specify our assumptions to the following points: 1) Solution for a complex arguments q 2 ,Imq 2 ^ is given 
by analytical continuation of Euclidean solution. 2) Minkowski solution corresponds with Imq 2 = 0, which is in 
principle obtainable by limiting procedure based on (1), can be directly obtained from continued Euclidean SDE, i.e. 
singularities are integrable and isolated. 3) The first singularity of the gluon propagator is located at the origin of 
complex plane. 

Continuation of PT-BFM gluon SDE is straightforward as there is only function g added when compared to the 
previous case. First let us rewrite Euclidean SDE by using Minkowski conventions q 2 E = —q 2 

Following the convention q 2 E — —q 2 for negative spacelike Minkowski q 2 = q 2 — q 2 then we can rewrite the Eq. (|3.4[) 
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which is valid for negative q 2 = q^ 
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i.e. now the integral domain on the rhs. of (|2.5|) becomes purely timelike as the all arguments appearing in SDE are. 
Again we have performed the substitution x — — z in order to show formal similarity with the Euclidean equation. 
The function g c in Eq. (|4.2I) should be the analytical continuation of our UV renormgroup improver, for which we 
take the function 



g 2 (z) = l/bln(e+\q 2 \/A 2 ) (4.4) 

where e = 2.73.. and we neglected in term as a legitimate simplification valid at ultraviolet. 

Here we briefly mention trivial fact that the Eq. (14.31) includes unphysical solutions, e.g. the real solution d(z) = 
dE(z), which corresponds to the choice <i -1 (0) = d^ 1 (0). We have also found some real solutions even for more 
"correct" choice d _1 (0) < 0. Again, they are clearly not hot candidates for wanted analytical continuation, as they 
do not agree with the dispersion relation (|2.3[) . 

We solve the Eq. (|4.2|) numerically at positive semiaxis q 2 — (0, +oo) by the method of iterations. We were 
searching for complex d which can have nonzero imaginary part everywhere at R + . Numerically we take e as constant 
satisfying e << |<i(0 + )| for all values of q 2 . Numerically e must be taken nonzero in order to accurately perform 
numerical integration in the vicinity of branch points in d. A bit unexpectedly a wide variety of stable numerical 
solutions of Eq (|4.2[) has been found. 

We anticipate here, the main feature of the solution is the appearance of two relatively large peaks with mutually 
opposite signs (see the results in figures 3-5) . Such singularities and the associated oscillations make the appropriate 
numerical findings a really hard problem. We achieve stable numerics by using Gaussian integrators such that intercept 
between two neighborhood points is smaller then e in the vicinity of singularities. We found advantageous to iterate 
not the function d, but rather its inverse which leads to very fast numerical convergence. Quite independently on the 
initial guess, typically ~ 30-40 iterations are enough to get vanishing (10~ 15 ) difference between last iterations. For 
instance the curves shown in Fig 2. were evaluated for 9000 Gaussian mesh points. 

In Euclidean space, the value of BFM gluon propagator at zero momenta represents a free parameter which must 
be fitted by hand. The value <i_E(0) = A was chosen for simplicity. There is no free choice when dealing with the SDE 
in Minkowski space and the right limit of d, ie. d(0 + ) value is necessarily constrained by global analytical property 
of d (stress here , we assume q 2 = is a branch point, in principle we allows the left and the right limit differ). Here 
we describe the numerical treatment which has been actually used in this paper. 

To find the best analytical continuation, we scan the values d(0 + ) and e, such that e/d(0 + ) << 1 and d(0 + ) ~ 
(0~) = A, for each set we find the solution of Minkowski SDEs. Then we check the dispersion relation for the solution 
Red(q 2 ) and Imd(q 2 ). The one which best reproduce the Euclidean solution via assumed dispersion relation for d 
is called the best analytical continuation (in sense of the point (2)). In accordance with spectral representation one 
expects vanishing imaginary part at the infrared as the correct analytical continuum limit. We can conclude, that 
within an accuracy limited by finite size of e one can always find the solution of timelike SDE. We expect the numerical 
accuracy of SDE solution is driven by the ratio, i.e e/A ~ 0.1 — 1%, while the agreement with assumed analyticity 
is a different issue. For this purpose we qualify the quality of analytical continuation by evaluating the following 
difference: 

-(q 2 ) = Q% (d E (q%) -I J <^^fp-) (4.5) 

where d E is the numerical solution of (|3.4j) and d is the numerical solution of (|4.3|) . The best analytical approximation 
was achieved by minimizing a and thus getting a few percentage deviation (for almost all q 2 ) at the minimum. For 
instance, for the case of small Schwinger coupling 7 — 0.1, the best analytical infrared fit we were able to trigger 
numerically uses e — 10~ 3 A 2 while keeping Red(Q + ) ~ d E (0) = A 2 . Let us stress the estimate of the numerical 
systematical error ~ 10 -3 (for absolute values of d) is always much smaller that the deviations from analyticity which 
can be estimated approximately by (d E (0) — d(0 + ))/ (d E (0) + d(0 + j) ~ 0.1 in the infrared. 



V. 7 DEPENDENCE OF MINKOWSKI SOLUTION 



The parameters 7 represents the strength of the effective coupling constant which arise due to the Schwinger 
mechanism. Thus numerical value of 7 should be responsible for the size of d(0) which can be estimated by the 
feedback of SDE solutions. Therefore we consider several values of 7 and report numerical results here for small 
7 = l/7r 2 and "large" 7 = ±1/tt couplings respectively. 
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In previous section we described the search of continued solution in details. In all studied cases we have found that 
the agreement with the exact analyticity is somehow limited ( a — in the exact case) . For a larger 7 it was difficult 
to get stable solution and minimize a so we were limited by concerning values | — ^ | < I. A few examples of a backward 
analytical continuation are shown in figure [1] for 7 = l/ir 2 and in figure [5] for 7 = — 1/tt, where they are compared 
with the original Euclidean solutions. Following the numerical procedure described, we label various solutions by 
the pair of numbers (Re c?~ 1 (0 + ), e) in the units of A. The difference between two initial choices (—0.902; 0.006) and 
(—0.9025; 0.006) exhibits numerical sensitivity in the case of small 7. 

While the dependence on 7 is rather mild in the case of the Euclidean solutions, we have found that this is not 
the case of the Minkowski space solutions. The solutions show up several cusps, number of them and their shapes 
largely change when varying 7. The only approximate invariant is the fact that the first maximum is positive (in 
sense of spectral function p) and located more or less at A, while the second peak is negative with the position shifted 
slightly to a higher scale. Other local maxima and minima are generated as well, also the discontinuity in the first 
derivatives are found. The first peak is located at scale A, which is in fact our choice of the infrared value of the 
gluon propagator, and quite expectingly, when the it is sharp enough, then one can see associated quasithreshold (at 
the scale 2A). The second peak is a surprise of the Minkowski solution. For a large momenta q 2 >> A 2 the second 
singularity regulates the contribution from the first one. Such screening is actually very large for small Schwinger 
coupling Considering the approximation p(uj) — S(uj — fi 2 ) — 5(ui — /1 2 ) of the infrared spectral function, the screening 
of the first peak by the second one is obvious. For instance for small 7 = 1/tt 2 we have found the position of the 
peaks are fii ~ A, \ii ~ V2.5A. 

In order to get some feeling how the numerical solution depends on the details which are not theoretically fully 
determined, we modify two parameters in the solely Minkowski SDE (|4.3[) , then after getting the solution we again 
use dispersion relation for d and see how the Euclidean solution is recovered. The main purpose to do this is to see 
whether the approximation employed has crucial influence on the analyticity. As a first, we have changed the sign of 7 
in the Minkowski SDE and then compare obtained solution with the original Euclidean one via backward continuation 
(|4.5|) . However alternating sign of 7 is ad hoc here, there exists a good motivation: beyond our truncation the constant 
7 will be replaced by the product of vertex function for glueball-like BSE for the gluon vertex. We expect , such 
product substantial changes when going to Minkowski timelike axis- it becomes complex and very likely enhanced 
there- therefore alternating the sign may not be so radical for our rough estimate. Actually, as we have found the 
deviations from analyticity a reaches smaller minima, such scenario is possible and perhaps preferred when compared 
to the case when 7 is unchanged. The resulting numerical solution are labeled as Model II in all presented figures, 
while what is labeled by letter I is calculated without any changes. Also the fits in figure 2 are based on the backward 
continuation of Minkowski solutions II which are displayed in figures 3-5. 

Secondly, we slightly change the improver g c , again solely in Minkowski SDE. For this purpose we replace the Euler 
constant e in the expression (I4.4|) by free parameter C. As a result, we have found g c enhanced (taking C = 2.5) is 
numerically preferred for positive value of 7, while taking g c slightly suppressed (taking C — 3.5) we can get better 
fits of Euclidean solution for negative 7. As an effect these changes modify the shape of the peaks. The appropriate 
solution for 7 = 1/tt coupling is added in figure [4] As one can see the second peak quite changes not only position 
but it transforms to a much narrower cusp. From this one can conclude that the analytical form of the resulting 
Minkowski propagator can be quite sensitive to the numerical approximations, which can be only marginal in the 
Euclidean treatment. 

As one can see from the Figs. 4 and 5, taking 7 small one get the solution with two very narrow peaks located 
at few A 2 . We could expect these peaks should be visible in a hadronic S- matrix element whenever timelike gluon is 
interchanged as an intermediate excitation. There would be "long living gluon" and "long living gluonic ghost" at the 
scale A according to their numerical "widths" which has the size by order smaller then A. Disregarding the reality 
of such scenario, we get interesting solution-two possible excitations appear for a single field A. Let us imagine what 
happens when we collide hadrons in such a world. Whenever energetically possible the quarks (and gluons itself) 
would preferably exchange low energy gluons with positive virtuality corresponding to the singularities or maxima 
positions and such peaks of gluonic colored corellator will become visible through the color subchannels. Even after 
the integration over the distribution of quarks, such gluonic colored mediators would be observed as resonances with 
very large cross sections (although without any doubt we know that many gluon interchanges are necessary to describe 
hadronic amplitudes, as only very large number can be responsible for formation of hadronic resonances, for very nice 
illustration see the study [38[ on tt — tt scattering). 

Increasing 7 changes the shape of the gluon propagator, especially the sharp peaks are melted down as they are 
transformed into more smooth maxima. For large 7 this is the branch point in the origin which is sole numerically 
observed singularity of the gluon propagator. The details of melting peaks is nicely seen in figure [SJ 
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VI. FITTING THE LATTICE 



The framework of PT-BFM SDEs is based on well denned theoretical background and in principle it could finally 
provide a gauge fixing independent set of new Greens functions. Nevertheless, the quantitative comparison with the 
standard-gauge fixed lattice simulation is further stimulating. Actually, there is small quantitative difference between 
decoupling Euclidean solutions, and the lattice evaluation of gluon propagator in Landau gauge [H, 5H, IHj] . In this 
respect, the analytical continuation obtained from PT-BFM SDE can be regarded as an reliable estimate of Minkowski 
space continuation of lattice data. 

Based on suitable numerical fits at Euclidean regime, a various form of continuation were proposed for scaling 
solutions of SDEs [44[ ■ According to the conclusion and numerical findings in [44| , the singularities of scaling solution 
studied are located near the real timelike semiaxis of the momenta. It is quite likely that the qualitative analytical 
properties of various solutions of YM SDEs systems are frankly insensitive to finiteness of d(0) and thus both, the 
scaling and the decoupling solutions of gluon propagators can be captured by the formula for Lehman representation. 
Then the quantitative difference will result from different forms of p, whatever it could be. Recent lattice data show 
finite infrared gluon propagator, which however is usual property of massive propagator, but as we have shown it 
does not necessary mean propagation of massive particle. There are legitimate attempts to move the first branch 
point towards the higher timelike scale [H], HI] , which fits are based on low q 2 lattice data. The first singularity fitted 
in this way assume quite obscured form of the spectral function. It is taken as a product of delta function and the 
inverse square root function singular at the point. Although it allows a reasonable estimate of the singularity position 
m = QOOMeV and it provides rough numerical fit of the Euclidean data at the region 0.5-5 GeV , the choice of the 
sing ularity form is far from obvious. The other lattice fits were proposed in the literature. For instance in the papers 



was used to describe the lattice prediction for the infrared gluon propagator up to k ~ l.5GeV. The function A has 
two complex conjugated poles and thus certainly does not belong to the sort of functions satisfying dispersion relation 
(|2.3j) . However, as it was shown recently in [40j, when A is considered in the loop, the resulting expression can be 
cast into the form of Eq. (|2.3|) (in fact, it is kind of an evidence, that albeit D does not have Lehman representation, 
its inverse has a part, which after the renormalization can be expressed through the dispersion relation). 

In our case, we have an additional hint estimated from the behavior of SDEs solution: the PT-BFM gluon propa- 
gator posses the first branch point at q = and the others are generated on the real axis by the Gauge Technique 
construction and selfconsistently determined by solving the continued SDE (|4.3[) . Making backward analytical contin- 
uation we reproduce the Euclidean solution approximately. From this we know, that we did not cross any singularities 
when one deforms usual Wick rotation contour at go plane. This is a strong argument which supports our findings, 
when comparing for instance to ad hoc guess of the singularity structure from fitting of shape of the Euclidean lattice 
data. 



The main purpose of our paper is to show that Stieltjes transformation or equivalently generalized Lehmann 
representation exists for PT-BFM gluon propagator. Very likely, the same applies for the conventional Landau gauge 
gluon propagator recently seen in the lattice simulation. We have found reliable spectral function for all q 2 , not 
only at the infrared or ultraviolet part separately. According to the expectations, the gluon propagator exhibits 
violation of positivity. In accordance with confinement the main singularity does not correspond with single delta 
function. However, for a small coupling characterizing the strength of the YM Schwinger-mechanism the infrared gluon 
propagator can be approximated by two mutually opposite sign Yukawa propagators, while the whole propagator 
function needs a certain amount of the continuous background in addition. Increasing 7 we get a more realistic 
picture of confinement -the sharp peaks in spectral function are melted down, while the Lehman approximation is 
less approximative, perhaps showing the evidence for complex conjugated poles of the from (|6.1J) . According to our 
approximative findings, the continuous spectrum starts at zero momenta, where we expect the first branch point. For 
small 7 , this is in excellent agreement with the assumption: the propagator is finite and purely real at light cone 
q 2 = 0. However for a larger 7 we get small imaginary part, which contradicts coexistence of Lehman representation 
and finiteness of the real part of gluon propagator at zero momenta. This, although small discrepancy tell us that 
our analytical assumption is not fulfilled , but is justified as an approximation only. 




A(k 2 ) 



k 2 + M 2 



(6.1) 



fc 4 + (M 2 + m 2 )k 2 + A 4 



VII. CONCLUSION 



0.1)1 



10(1 



L IK* 

q E /A 



FIG. 1: Spacelike solution for the model with coupling 7 = 1/iv 2 and d _1 (CP) = A 2 . Comparison with spectral representation 
fits are added, the bracket shows the numerical fit: (i?ed(0 + )~ 1 , e) (in units of A 2 ) for which the Minkowski SDE has been 
solved. 



In order to discuss a possible sources of this discrepance let us tentatively divide the gluon PT-BFM propagator 
into two peaces: 



where tna{<1 2 ) should be small when compared to the full d{q 2 ) and we distinguish among several possibilities. 

1. where r^A is the marginal complex remnant which does not obey Stieltjes representation but does not contradict 
Wick rotation and thus the primacy of Euclidean calculation as well. The example is complex pole located at point 
z with Rex > and Imz < 

2. rpjA does not allow Wick rotation, there are singularities and cuts in complex plane of q 2 (and qo) which makes 
the theories in Euclidean and Minkowski space different. In this case, the lattice data should be regarded as a certain 
approximation of Minkowski spacelike correlators. The example of such behaviour we can mention is the function 
given by (|6.1[) , if it has complex pole at some z, Rez > 0, Imz > 0. 

3- rpfA vanishes for the full solution. It is possible it is a consequence of approximation employed, eg. it is a possible 
artifact of weakness of the Gauge Technique linearization, truncation of SDE system, modeling Schwinger mechanism 
and modeling UV RG improver g c . In this case the idea of analytical effective charge finds new applications beyond 
its original perturbative conjecture [45l - |48j 

To distinguish among above three points is beyond our recent analyzing power, however in future study one should 
always ask, what improvement of the knowledge has been achieved in this respect. 

In practice, most of theoretical description of QCD processes is based on the factorization schemes where the hard 
part of amplitudes is calculated by using perturbation theory while the non-perturbative information is included in 
the partonic distribution functions. Therefore the first hint we can get, can be obtained from the phenomenological 
studies which adopt the running charges constrained by finite BFM gluon propagator, for recent calculation of the 
proton structure function see (3lj . and [32[ for heavy meson decay study. The timelike gluons exchanged certainly 
appear in the s-channel of meson scattering or heavy meson strong decays. The calculation difficulty of such process is 
obvious, the observed gluon peaks will interfere not only with itself, but with the first colorless resonances, e.g. p meson 
which always arise nonperturbatively when very large number (infinite) gluons are exchanged ( 38]). Nevertheless the 
difficulty of the problem, we believe that propagator calculated here can be useful when one would like to investigate 
glueball and hadron spectra in cases when timelike kinematics is especially pronounced. For this purpose inclusion of 




(7.1) 
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A f=lk 

fit J=V% ; A"'(OH-0.66;0.063) 
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0.01 



ion 
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FIG. 2: Spacelike solution for various 7. Dashed and dotted lines represent two fits obtained for —7 = 1/tt from Minkowski 
SDE solution (see the text for the details). The fits are labeled by values of Red(0 + )~ 1 , e in units of A 2 . 



the quark fields remains to be done. 
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FIG. 4: Magnitude of PT-BFM gluon propagator in Minkowski space for Model I: 7 = 1/tt 2 and 7?ed(0 + ) _1 = 0.9A, e = 0.006A 2 , 
corresponding d^J 1 (0) = A 2 and the Model II, which is the same as in the Fig. [3] Dotted lines represents Model II but with 
the function g c as described in the text. 




FIG. 5: Details of the singularities (and shape) of the timelike gluon form factor as obtained for small 7 = 1/-7T 2 (Model I) and 
moderate coupling 7 = 1/tt (Model II) respectively. 



